Correction: Entamoeba lysyl-tRNA Synthetase Contains a Cytokine-Like Domain with Chemokine Activity towards Human Endothelial Cells.

[This corrects the article DOI: 10.1371/journal.pntd.0001398.].

Cellular and clinical impact of protein phosphatase enzyme epigenetic silencing in multiple cancer tissues.

BACKGROUND: Protein Phosphatase Enzymes (PPE) and protein kinases simultaneously control phosphorylation mechanisms that tightly regulate intracellular signalling pathways and stimulate cellular responses. In human malignancies, PPE and protein kinases are frequently mutated resulting in uncontrolled kinase activity and PPE suppression, leading to cell proliferation, migration and resistance to anti-cancer therapies. Cancer associated DNA hypermethylation at PPE promoters gives rise to transcriptional silencing (epimutations) and is a hallmark of cancer. Despite recent advances in sequencing technologies, data availability and computational capabilities, only a fraction of PPE have been reported as transcriptionally inactive as a consequence of epimutations. METHODS: In this study, we examined promoter-associated DNA methylation profiles in Protein Phosphatase Enzymes and their Interacting Proteins (PPEIP) in a cohort of 705 cancer patients in five tissues (Large intestine, Oesophagus, Lung, Pancreas and Stomach) in three cell models (primary tumours, cancer cell lines and 3D embedded cancer cell cultures). As a subset of PPEIP are known tumour suppressor genes, we analysed the impact of PPEIP promoter hypermethylation marks on gene expression, cellular networks and in a clinical setting. RESULTS: Here, we report epimutations in PPEIP are a frequent occurrence in the cancer genome and manifest independent of transcriptional activity. We observed that different tumours have varying susceptibility to epimutations and identify specific cellular signalling networks that are primarily affected by epimutations. Additionally, RNA-seq analysis showed the negative impact of epimutations on most (not all) Protein Tyrosine Phosphatase transcription. Finally, we detected novel clinical biomarkers that inform on patient mortality and anti-cancer treatment sensitivity. CONCLUSIONS: We propose that DNA hypermethylation marks at PPEIP frequently contribute to the pathogenesis of malignancies and within the precision medicine space, hold promise as biomarkers to inform on clinical features such as patient survival and therapeutic response.

ERBB4-mediated signaling is a mediator of resistance to BTK and PI3K inhibitors in B cell lymphoid neoplasms.

BTK and PI3K inhibitors are among the drugs approved for the treatment of patients with lymphoid neoplasms. Although active, their ability to lead as single agents to long-lasting complete remission is rather limited especially in the lymphoma setting. This indicates that tumor cells often develop resistance to the drugs. Here, we show that the overexpression of ERBB4 and its ligands represents a modality for B cell neoplastic cells to bypass the anti-tumor activity of BTK and PI3K inhibitors and that targeted pharmacological interventions can restore sensitivity to the small molecules. We started from a marginal zone lymphoma (MZL) cell line, Karpas-1718, kept under prolonged exposure to the PI3Kδ inhibitor idelalisib until acquisition of resistance, or with no drug. Cells underwent transcriptome, miRNA and methylation profiling, whole exome sequencing, and pharmacological screening which led to the identification of the overexpression of ERBB4 and its ligands HBEGF and NRG2 in the resistant cells. Cellular and genetic experiments demonstrated the involvement of this axis in blocking the anti-tumor activity of various BTK and PI3K inhibitors, currently used in the clinical setting. Addition of recombinant HBEGF induced resistance to BTK and PI3K inhibitors in parental cells but also in additional lymphoma models. Combination with the ERBB inhibitor lapatinib was beneficial in resistant cells and in other lymphoma models already expressing the identified resistance factors. Multi-omics analysis underlined that an epigenetic reprogramming affected the expression of the resistance-related factors, and pretreatment with demethylating agents or EZH2 inhibitors overcame the resistance. Resistance factors were shown to be expressed in clinical samples, further extending the findings of the study. In conclusions, we identified a novel ERBB4-driven mechanism of resistance to BTK and PI3K inhibitors and treatments that appear to overcome it. Key points: A mechanism of secondary resistance to the PI3Kδ and BTK inhibitors in B cell neoplasms driven by secreted factors.Resistance can be reverted by targeting ERBB signaling.

Genetically predicted telomere length and Alzheimer's disease endophenotypes: a Mendelian randomization study.

Telomere length (TL) is associated with biological aging, consequently influencing the risk of age-related diseases such as Alzheimer's disease (AD). We aimed to evaluate the potential causal role of TL in AD endophenotypes (i.e., cognitive performance, N = 2233; brain age and AD-related signatures, N = 1134; and cerebrospinal fluid biomarkers (CSF) of AD and neurodegeneration, N = 304) through a Mendelian randomization (MR) analysis. Our analysis was conducted in the context of the ALFA (ALzheimer and FAmilies) study, a population of cognitively healthy individuals at risk of AD. A total of 20 single nucleotide polymorphisms associated with TL were used to determine the effect of TL on AD endophenotypes. Analyses were adjusted by age, sex, and years of education. Stratified analyses by APOE-ɛ4 status and polygenic risk score of AD were conducted. MR analysis revealed significant associations between genetically predicted longer TL and lower levels of CSF Aß and higher levels of CSF NfL only in APOE-ɛ4 non-carriers. Moreover, inheriting longer TL was associated with greater cortical thickness in age and AD-related brain signatures and lower levels of CSF p-tau among individuals at a high genetic predisposition to AD. Further observational analyses are warranted to better understand these associations.

Altered methylation pattern in EXOC4 is associated with stroke outcome: an epigenome-wide association study.

BACKGROUND AND PURPOSE: The neurological course after stroke is highly variable and is determined by demographic, clinical and genetic factors. However, other heritable factors such as epigenetic DNA methylation could play a role in neurological changes after stroke. METHODS: We performed a three-stage epigenome-wide association study to evaluate DNA methylation associated with the difference between the National Institutes of Health Stroke Scale (NIHSS) at baseline and at discharge (ΔNIHSS) in ischaemic stroke patients. DNA methylation data in the Discovery (n = 643) and Replication (n = 62) Cohorts were interrogated with the 450 K and EPIC BeadChip. Nominal CpG sites from the Discovery (p value < 10-06) were also evaluated in a meta-analysis of the Discovery and Replication cohorts, using a random-fixed effect model. Metabolic pathway enrichment was calculated with methylGSA. We integrated the methylation data with 1305 plasma protein expression levels measured by SOMAscan in 46 subjects and measured RNA expression with RT-PCR in a subgroup of 13 subjects. Specific cell-type methylation was assessed using EpiDISH. RESULTS: The meta-analysis revealed an epigenome-wide significant association in EXOC4 (p value = 8.4 × 10-08) and in MERTK (p value = 1.56 × 10-07). Only the methylation in EXOC4 was also associated in the Discovery and in the Replication Cohorts (p value = 1.14 × 10-06 and p value = 1.3 × 10-02, respectively). EXOC4 methylation negatively correlated with the long-term outcome (coefficient = - 4.91) and showed a tendency towards a decrease in EXOC4 expression (rho = - 0.469, p value = 0.091). Pathway enrichment from the meta-analysis revealed significant associations related to the endocytosis and deubiquitination processes. Seventy-nine plasma proteins were differentially expressed in association with EXOC4 methylation. Pathway analysis of these proteins showed an enrichment in natural killer (NK) cell activation. The cell-type methylation analysis in blood also revealed a differential methylation in NK cells. CONCLUSIONS: DNA methylation of EXOC4 is associated with a worse neurological course after stroke. The results indicate a potential modulation of pathways involving endocytosis and NK cells regulation.

DDR1 methylation is associated with bipolar disorder and the isoform expression and methylation of myelin genes.

Aim: To investigate DDR1 methylation in the brains of bipolar disorder (BD) patients and its association with DDR1 mRNA levels and comethylation with myelin genes. Materials & methods: Genome-wide profiling of DNA methylation (Infinium MethylationEPIC BeadChip) corrected for glial composition and DDR1 gene expression analysis in the occipital cortices of individuals with BD (n = 15) and healthy controls (n = 15) were conducted. Results:DDR1 5-methylcytosine levels were increased and directly associated with DDR1b mRNA expression in the brains of BD patients. We also observed that DDR1 was comethylated with a group of myelin genes. Conclusion:DDR1 is hypermethylated in BD brain tissue and is associated with isoform expression. Additionally, DDR1 comethylation with myelin genes supports the role of this receptor in myelination.

Publisher Correction: MethCORR modelling of methylomes from formalin-fixed paraffin-embedded tissue enables characterization and prognostication of colorectal cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

DNA methylation of MMPs and TIMPs in atherothrombosis process in carotid plaques and blood tissues.

BACKGROUND AND PURPOSE: Polymorphisms and serum levels of Matrix Metalloproteinases (MMP) and Tissue Inhibitor of Metalloproteinases (TIMP) have been studied with regard to atheromatous plaques and ischemic stroke, while no studies of DNA methylation (DNAm) patterns of MMP or TIMP have been performed to that end. Here, we evaluate DNAm levels of the MMP and TIMP gene families in human carotid plaques and blood samples of atherothrombotic stroke patients. METHODS: We profiled the DNAm status of stable and ulcerated atherosclerotic plaques obtained as pair sets from three patients who underwent carotid endarterectomy surgery. We selected 415 CpG sites, mapping into MMPs and TIMPs genes for further study. Secondly, the statistically associated CpG sites were analyzed in blood samples from two separate atherothrombotic stroke cohorts (total sample size = 307), ischemic stroke-cohort 1 (ISC-1): 37 atherothrombotic patients and 6 controls, ischemic stroke-cohort 2 (ISC-2): 80 atherothrombotic patients and 184 controls. DNAm levels from plaque tissue and blood samples were evaluated using a high-density microarray Infinium, HumanMethylation450 BeadChip and Infinium MethylationEPIC BeadChip. RESULTS: Three CpG sites were statistically significantly associated with unstable plaque portions; cg02969624, q-value = 0.035 (TIMP2), and cg04316754, q-value = 0.037 (MMP24) were hypermethylated, while cg24211657 q-value = 0.035 (TIMP2) was hypomethylated. Association of cg04316754 (MMP24) methylation levels with atherothrombotic risk was also observed in blood tissue: ISC-1 p-values = 0.03, ISC-2 p-value = 1.9 × 10-04. CONCLUSIONS: The results suggest different DNAm status of MMP24 between stable and unstable atherothrombotic carotid plaques, and between atherothrombotic stroke and controls in blood samples.

Epigenetic loss of RNA-methyltransferase NSUN5 in glioma targets ribosomes to drive a stress adaptive translational program.

Tumors have aberrant proteomes that often do not match their corresponding transcriptome profiles. One possible cause of this discrepancy is the existence of aberrant RNA modification landscapes in the so-called epitranscriptome. Here, we report that human glioma cells undergo DNA methylation-associated epigenetic silencing of NSUN5, a candidate RNA methyltransferase for 5-methylcytosine. In this setting, NSUN5 exhibits tumor-suppressor characteristics in vivo glioma models. We also found that NSUN5 loss generates an unmethylated status at the C3782 position of 28S rRNA that drives an overall depletion of protein synthesis, and leads to the emergence of an adaptive translational program for survival under conditions of cellular stress. Interestingly, NSUN5 epigenetic inactivation also renders these gliomas sensitive to bioactivatable substrates of the stress-related enzyme NQO1. Most importantly, NSUN5 epigenetic inactivation is a hallmark of glioma patients with long-term survival for this otherwise devastating disease.

Epigenetic loss of the endoplasmic reticulum-associated degradation inhibitor SVIP induces cancer cell metabolic reprogramming.

The endoplasmic reticulum (ER) of cancer cells needs to adapt to the enhanced proteotoxic stress associated with the accumulation of unfolded, misfolded and transformation-associated proteins. One way by which tumors thrive in the context of ER stress is by promoting ER-Associated Degradation (ERAD), although the mechanisms are poorly understood. Here, we show that the Small p97/VCP Interacting Protein (SVIP), an endogenous inhibitor of ERAD, undergoes DNA hypermethylation-associated silencing in tumorigenesis to achieve this goal. SVIP exhibits tumor suppressor features and its recovery is associated with increased ER stress and growth inhibition. Proteomic and metabolomic analyses show that cancer cells with epigenetic loss of SVIP are depleted in mitochondrial enzymes and oxidative respiration activity. This phenotype is reverted upon SVIP restoration. The dependence of SVIP hypermethylated cancer cells on aerobic glycolysis and glucose was also associated with sensitivity to an inhibitor of the glucose transporter GLUT1. This could be relevant to the management of tumors carrying SVIP epigenetic loss, because these occur in high-risk patients who manifest poor clinical outcomes. Overall, our study provides insights into how epigenetics helps deal with ER stress and how SVIP epigenetic loss in cancer may be amenable to therapies that target glucose transporters.

Circular RNA CpG island hypermethylation-associated silencing in human cancer.

Noncoding RNAs (ncRNAs), such as microRNAs and long noncoding RNAs (lncRNAs), participate in cellular transformation. Work done in the last decade has also demonstrated that ncRNAs with growth-inhibitory functions can undergo promoter CpG island hypermethylation-associated silencing in tumorigenesis. Herein, we wondered whether circular RNAs (circRNAs), a type of RNA transcripts lacking 5'-3' ends and forming closed loops that are gaining relevance in cancer biology, are also a target of epigenetic inactivation in tumors. To tackle this issue, we have used cancer cells genetically deficient for the DNA methyltransferase enzymes in conjuction with circRNA expression microarrays. We have found that the loss of DNA methylation provokes a release of circRNA silencing. In particular, we have identified that promoter CpG island hypermethylation of the genes TUSC3 (tumor suppressor candidate 3), POMT1 (protein O-mannosyltransferase 1), ATRNL1 (attractin-like 1) and SAMD4A (sterile alpha motif domain containing 4A) is linked to the transcriptional downregulation of both linear mRNA and the hosted circRNA. Although some circRNAs regulate the linear transcript, we did not observe changes in TUSC3 mRNA levels upon TUSC3 circ104557 overexpression. Interestingly, we found circRNA-mediated regulation of target miRNAs and an in vivo growth inhibitory effect upon TUSC3 circ104557 transduction. Data mining for 5'-end CpG island methylation of TUSC3, ATRNL1, POMT1 and SAMD4A in cancer cell lines and primary tumors showed that the epigenetic defect was commonly observed among different tumor types in association with the diminished expression of the corresponding transcript. Our findings support a role for circRNA DNA methylation-associated loss in human cancer.

Bromodomain inhibition shows antitumoral activity in mice and human luminal breast cancer.

BET bromodomain inhibitors, which have an antitumoral effect against various solid cancer tumor types, have not been studied in detail in luminal breast cancer, despite the prevalence of this subtype of mammary malignancy. Here we demonstrate that the BET bromodomain inhibitor JQ1 exerts growth-inhibitory activity in human luminal breast cancer cell lines associated with a depletion of the C-MYC oncogene, but does not alter the expression levels of the BRD4 bromodomain protein. Interestingly, expression microarray analyses indicate that, upon JQ1 administration, the antitumoral phenotype also involves downregulation of relevant breast cancer oncogenes such as the Breast Carcinoma-Amplified Sequence 1 (BCAS1) and the PDZ Domain-Containing 1 (PDZK1). We have also applied these in vitro findings in an in vivo model by studying a transgenic mouse model representing the luminal B subtype of breast cancer, the MMTV-PyMT, in which the mouse mammary tumor virus promoter is used to drive the expression of the polyoma virus middle T-antigen to the mammary gland. We have observed that the use of the BET bromodomain inhibitor for the treatment of established breast neoplasms developed in the MMTV-PyMT model shows antitumor potential. Most importantly, if JQ1 is given before the expected time of tumor detection in the MMTV-PyMT mice, it retards the onset of the disease and increases the survival of these animals. Thus, our findings indicate that the use of bromodomain inhibitors is of great potential in the treatment of luminal breast cancer and merits further investigation.

Epigenetic profiling to classify cancer of unknown primary: a multicentre, retrospective analysis.

BACKGROUND: Cancer of unknown primary ranks in the top ten cancer presentations and has an extremely poor prognosis. Identification of the primary tumour and development of a tailored site-specific therapy could improve the survival of these patients. We examined the feasability of using DNA methylation profiles to determine the occult original cancer in cases of cancer of unknown primary. METHODS: We established a classifier of cancer type based on the microarray DNA methylation signatures (EPICUP) in a training set of 2790 tumour samples of known origin representing 38 tumour types and including 85 metastases. To validate the classifier, we used an independent set of 7691 known tumour samples from the same tumour types that included 534 metastases. We applied the developed diagnostic test to predict the tumour type of 216 well-characterised cases of cancer of unknown primary. We validated the accuracy of the predictions from the EPICUP assay using autopsy examination, follow-up for subsequent clinical detection of the primary sites months after the initial presentation, light microscopy, and comprehensive immunohistochemistry profiling. FINDINGS: The tumour type classifier based on the DNA methylation profiles showed a 99·6% specificity (95% CI 99·5-99·7), 97·7% sensitivity (96·1-99·2), 88·6% positive predictive value (85·8-91·3), and 99·9% negative predictive value (99·9-100·0) in the validation set of 7691 tumours. DNA methylation profiling predicted a primary cancer of origin in 188 (87%) of 216 patients with cancer with unknown primary. Patients with EPICUP diagnoses who received a tumour type-specific therapy showed improved overall survival compared with that in patients who received empiric therapy (hazard ratio [HR] 3·24, p=0·0051 [95% CI 1·42-7·38]; log-rank p=0·0029). INTERPRETATION: We show that the development of a DNA methylation based assay can significantly improve diagnoses of cancer of unknown primary and guide more precise therapies associated with better outcomes. Epigenetic profiling could be a useful approach to unmask the original primary tumour site of cancer of unknown primary cases and a step towards the improvement of the clinical management of these patients. FUNDING: European Research Council (ERC), Cellex Foundation, the Institute of Health Carlos III (ISCIII), Cancer Australia, Victorian Cancer Agency, Samuel Waxman Cancer Research Foundation, the Health and Science Departments of the Generalitat de Catalunya, and Ferrer.

Analysis of gene expression changes in peach leaves in response to Plum pox virus infection using RNA-Seq.

Differences in gene expression were studied after Plum pox virus (PPV, sharka disease) infection in peach GF305 leaves with and without sharka symptoms using RNA-Seq. For each sample, more than 80% of 100-nucleotide paired-end (PE) Illumina reads were aligned on the peach reference genome. In the symptomatic sample, a significant proportion of reads were mapped to PPV reference genomes (1.04% compared with 0.00002% in non-symptomatic leaves), allowing for the ultra-deep assembly of the complete genome of the PPV isolate used (9775 nucleotides, missing only 11 nucleotides at the 5' genome end). In addition, significant alternative splicing events were detected in 359 genes and 12 990 single nucleotide polymorphisms (SNPs) were identified, 425 of which could be annotated. Gene ontology annotation revealed that the high-ranking mRNA target genes associated with the expression of sharka symptoms are mainly related to the response to biotic stimuli, to lipid and carbohydrate metabolism and to the negative regulation of catalytic activity. A greater number of differentially expressed genes were observed in the early asymptomatic phase of PPV infection in comparison with the symptomatic phase. These early infection events were associated with the induction of genes related to pathogen resistance, such as jasmonic acid, chitinases, cytokinin glucosyl transferases and Lys-M proteins. Once the virus had accumulated, the overexpression of Dicer protein 2a genes suggested a gene silencing plant response that was suppressed by the virus HCPro and P1 proteins. These results illustrate the dynamic nature of the peach-PPV interaction at the transcriptome level and confirm that sharka symptom expression is a complex process that can be understood on the basis of changes in plant gene expression.